Vortex lines or sheets — what is formed in dynamic drives? 
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In isotropic macroscopic quantum systems vortex lines can be formed while in anisotropic systems 
also vortex sheets are possible. Based on measurements of superfluid ''He-A, we present the principles 
which select between these two competing forms of quantized vorticity: sheets displace lines if the 
frequency of the external fleld exceeds a critical limit. The resulting topologically stable state 
consists of multiple vortex sheets and has much faster dynamics than the state with vortex lines. 



The state of superfluid ^He in rotation was originally 
explained by Onsager and Landau and Lifshitz by 
concentrating the required velocity circulation into co- 
axial sheets which, when viewed on large scales, give rise 
to solid-body-like rotation. Later, Feynman |^ showed 
that sheets are not stable in an isotropic superfluid and 
that quantized vortex lines are formed instead. The same 
explanation was provided by Abrikosov for type II 
superconductors in a magnetic field. 

Recently, however, the vortex sheet has staged a come 
back: It was theoretically suggested for unconventional 
superconductors |5j and experimentally verified in the 
anisotropic p-wave superfluid '^He-A The sheets are 
formed from a topologically stable domain-wall-like pla- 
nar defect into which the vorticity is confined. Since these 
quantum systems also support vortex lines, the contro- 
versy remains: What determines the formation of lines 
or sheets, or alternatively, what is the difference in the 
macroscopic physical properties of the system when there 
are lines or sheets? This question is not limited to super- 
conductors or '^He-A, but is of importance for defect for- 
mation in general, possibly in the early universe, in neu- 
tron stars, Bose-Einstein condensates, etc. - wherever 
non-trivial symmetry breaking into a multi-component 
order-parameter field may occur. 

The answer whether lines or sheets are formed de- 
pends, according to our measurements of '^He-A, on the 
applied field: If its frequency and amplitude exceed the 
appropriate critical limits, sheets appear and displace 
lines, as shown schematically in Fig. ^ The measured 
phase diagram is supported by our analysis of the dy- 
namic properties of vortex sheets, and is further con- 
firmed by numerical simulations. The transition from 
lines to sheets is driven by the much faster dynamics of 
states with multiple vortex sheets. The roots of this phe- 
nomenon share common features with self organization 
and pattern formation in classical systems. 

Vortex lines vs sheets.- Superconductors respond to 



FIG. 1. The different topologically distinct states of ro- 
tating ^He-A depend on the amplitude and frequency of the 
drive: (a) The vortex-free Meissner state exists below a crit- 
ical value of the amplitude, (b) The state with quantized 
vortex lines is created when the drive increases sufficiently 
slowly, (c) Vortex sheets are created at higher frequencies. 
The critical values depend on the global order-parameter tex- 
ture and display large metastability. 



a magnetic field by producing screening currents. Sim- 
ilarly, a superfluid attempts to corotate at the con- 
tainer's angular velocity fl. The two motions resem- 
ble each other and usually take the form of quantized 
vortex lines so that the superfluid vorticity V x is 
non-zero only within vortex cores, having on average the 
solid-body value 2fl. In a state with sheets the irrota- 
tional superflow is channelled parallel to a sheet, while 
Vs jumps across the sheet leaving the non-zero vortic- 
ity within the sheet. The distance between the sheets is 

1 /3 

b — (Sct/ps^^^) , where a is the surface tension of the 
sheet and Ps the superfluid density [^j. In equilibrium, 
vortex lines in '^He-A have a lower energy than vortex 
sheets 0. 

The stability of the vortex sheet rests on the existence 
of a two-fold degeneracy of the ground state, so that a 
topologically stable domain wall can exist. In '^He-A, 
it separates energetically equivalent states with parallel 
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and antiparallel orientations of the orbital 1 and spin d 
anisotropy directions. In superconductors with broken 
time-reversal symmetry the degeneracy could arise from 
the locking of 1 to the crystal- field anisotropy |5j. The 
vorticity in the domain wall is quantized it is dis- 
tributed along the wall with a period p — nn / {2U,h) . Here 
K is the circulation quantum and n the topological charge 
which depends on the structure of the vorticity. For the 
vortex sheet in ^He-A n — \ while separate vortex lines 
are generally doubly quantized {n — 2) 

The large-scale structure of the vortex sheet may take 
topologically distinct configurations with different num- 
bers M of separated sheets. The minimum is one contin- 
uous sheet {M = 1) in the form of a double spiral cen- 
tered around the container axis (see bottom-right inset 
in Fig. H) ; this is the lowest energy state of a vortex sheet 
at constant fJ. This configuration is closest to the con- 
centric vortex sheets originally envisioned by Onsager 0] 
and by Landau and Lifshitz ||] . The other extreme is the 
maximum number of sheets, Mmax = 7ri?/&, which corre- 
sponds to one connection with the lateral sample bound- 
ary per each inter-layer distance b along the perimeter 
of the sample (Mmax « 20 at = Irad/s). In the ab- 
sence of pinning, like in ■^He-A, the sheets fold to fill the 
sample evenly since a uniform distribution of vorticity 
is energetically preferred. Mutual repulsion between the 
sheets and the boundary conditions force the sheets to be 
oriented perpendicular to the boundary. This promotes 
their radial orientation, as depicted in Fig. ^ (c). 

Dynamics of lines and sheets - A time-dependent ro- 
tation drive 17 (i) exerts a force on the circulation quanta 
which is mostly radial in solid-body-like rotation. The 
vorticity averaged over the sample, fii, — (|V x Vs|)/2, 
obeys (§ 

da„/dt + 2aVL^ {VL^ - f^) = 0, (1) 

where a is the mobility in the radial direction determined 
by the friction from the normal excitations. Owing to 
the strongly anisotropic order-parameter orientations in 
the sheet, the mobility of vorticity in the sheet is also 
anisotropic. It is determined by the spatial inhomogene- 
ity of the orbital anisotropy vector \ : 

ar'^a,' J d\{v,l)' . (2) 

Here i refers to the directions along || or perpendicular 
_L to the sheet, and ao = pnPsi^/ i'^P/J-) is related to the 
Cross- Anderson orbital viscosity p, We neglect an 

order-of-one anisotropy of ps- One obtains ay ~ ao{p/s) 
and a± ^ uq^s/p), where s w 40 /im is the thickness 
of the sheet. The mobility of separate vortex lines is 
ay ~ ckQ , since the dimensions of the vortex core are s 
in all directions. Thus the mobility along the sheet 
is the highest: a|| > > a±. At 17 = 1 rad/s one has 
a II : a-u : a± ~ 4 : 1 : 1/4. Since it is the radial mobility 



that matters in Eq. (|l|), radial orientation of the sheets 
improves their dynamic response. 

For the dynamics also the creation and annihilation of 
vorticity is crucial. The circulation quanta are created at 
a critical velocity which depends on the order-parameter 
texture ||l2|. For vortex sheets, vorticity is created at the 
connection lines of the sheets with the boundary where 
the order-parameter field is already distorted and the 
critical velocity is reduced compared to vortex lines 
Thus multiple sheets with their 2AI connection lines pro- 
vide a large number of effective nucleation sites. These 
inject new vorticity much more efficiently than the com- 
peting alternative, one nucleation center of vortex lines 
close to the lateral sample boundary 

To conclude, the different rates of dynamic response, 
which control relaxation towards the instantaneous equi- 
librium, are expected to govern the transition between 
lines and sheets: A high-frequency drive at sufficient am- 
plitude should replace all other forms of vorticity with ra- 
dially oriented multiple vortex sheets. The transition can 
occur in the following way: At high frequency the order- 
parameter texture around a nucleation center, where vor- 
tex lines are periodically nucleated and annihilated, may 
become locally unstable and the beginning of a domain- 
wall defect may be created. When existing vortex lines 
have been annihilated during a subsequent decrease in 
the drive, the new vorticity may now re-enter during the 
next increase in the drive via the domain- wall defect, 
where the critical velocity is reduced and the mobility 
is enhanced. A single sheet is unstable in an oscillating 
drive and, as discussed later, this leads to the formation 
of a state with multiple sheets. A configuration with a 
given number of sheets is topologically stable and will 
persist as a metastable state if the time-dependent com- 
ponent in the drive is switched off and Q is kept constant. 

Measurement of topological transitions - Our measure- 
ments on ^He-A have been performed on a cylindrical 
sample with a radius R — 2 mm in an axial magnetic 
field of 9.6 mT, pressure 33 bar, and temperature 1.7 mK, 
by rotating the entire refrigerator. NMR spectroscopy 
is used to identify the various order-parameter topolo- 
gies. Vortex lines or sheets produce characteristic NMR 
satellite peaks, while in the vortex-free state there is no 
satellite. Although the satellites may overlap (Fig. 2, in- 
sert), they have different dependencies on il which distin- 
guishes between lines and sheets |6|. The rotation drive 
is thus used both for generating transitions between the 
states and as a detection tool. 

The lower part of Fig. ^ illustrates transitions between 
the three rotating states of Fig. |l| as a function of time 
t when the drive Q{t) ~ ^\smujt is switched on. The 
starting state is vortex free with no satellite. Vortex lines 
begin to form when 17 reaches the critical value f7c ~ 
0.24 rad/s during the first half period. Here the height of 
the satellite is proportional to the total number of vortex 
lines. The characteristic features are two plateaus during 
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FIG. 2. (Top) Sequence of transitions (vortex-free vor- 
tex lines — > vortex sheets), when the drive Q{t) = Qisinut 
(with Qi = 0.35 rad/s) is switched on: (•) first appearance of 
lines, (0) transition to sheets starts, (a) transition to sheets is 
completed. The vertical axis is the period 27r/a; of the drive 
while the horizontal axis is the common time axis (t) with 
the bottom panel. (Bottom) Height of the vorticity satellites 
as a function of time after switching on the drive fl(t), with 
2n/Lj = 600 s which corresponds to the path along the dashed 
arrow in the top panel. The signal is recorded at the NMR 
frequency shown in the insert. The drive fi(t) is plotted below 
the satellite signal. (Insert) NMR absorption spectra of the 
the rotating states with lines and sheets at Q — 2.4 rad/s and 
of the vortex- free state (Q < fie). 

one half cycle: the first at zero height when \Q\ increases 
from to Qc, and the second at maximum absorption 
vifhcn decreases from fti to roughly fli — ilc- The 
plateaus correspond to zero and to maximum number of 
vortex lines, respectively. 

After one full period the conversion to vortex sheets 
starts. The response gradually changes and acquires fea- 
tures characteristic to the vortex sheet: The satellite 
height grows larger and the plateaus disappear. For a 
sheet, the signal height is proportional to its total length. 
The critical velocity is now vanishingly small which en- 
sures the smooth growth and shrinkage of the sheets as 
a function of fl. Also, the satellite height does not reach 
zero when fl changes sign because the backbone of the 
sheet - the domain wall - does not have time to drift 
away and continues to contribute to the signal even in 
the absence of circulation at SI = 0. These features can 
be studied in detail, by switching off the modulation of 
the drive at any moment. Thereafter fl can be kept at 
constant value, and the entire NMR absorption spectrum 
can be recorded to determine the frequency shift and 
height of the satellite. 

The upper left panel of Fig. || illustrates how long it 
takes to complete the transitions at different u. Note 
that there exists a characteristic frequency uj* such that 
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FIG. 3. Measured hysteresis loops in satellite-peak height 
when Q is slowly changed (solid lines), (a) An adiabati- 
cally grown single sheet, (b) Multiple sheets created from 
the single sheet by subjecting it to high-frequency modula- 
tion [n{t) = (0.4 smut+ 1.2) rad/s, with 2tt/lj = 10s]. The 
broken lines represent numerical simulations for the config- 
urations shown in the respective inserts. The only fitting 
parameter is the coefficient of proportionality which connects 
the total length of the sheets to NMR absorption. 

2tt/uj* ~ 10'^ s which separates lines and sheets in the 
limit t —^ oo: At to < uj* a response in terms of lines is 
stable while at w > w* sheets appear. The conditions 
in the lower panel correspond to the case when u; > uj* 
and sheets are stable. The larger the drive frequency lu, 
the faster the transition to sheets starts. Finally, with 
2tt/uj < 20 s the transition to sheets starts directly from 
the vortex-free state, with no indication of a signal from 
vortex lines in between. With increasing drive amplitude 
rii the transition to sheets oj* moves to lower frequencies 
and is completed faster. Thus the transition line u!*{il.i) 
is tilted as shown in Fig. |l|. 

Identification of sheet topologies.- Different vortex- 
sheet configurations can be studied by measuring the f2- 
dependence of the NMR satellite. This can be recorded 
either in the adiabatically slow limit (Fig. ^) or as a dy- 
namic response (Fig. To analyze the response of a 
state with multiple sheets it is instructive to compare it 
to that of the well-known configuration with a single con- 
tinuously coiled sheet. The latter provides a convenient 
reference which can be grown adiabatically from a single 
seed of domain wall, by slowly increasing from zero 

Let us calculate the total sheet length L for a given 
number of sheets M and circulation quanta N. The vor- 
ticity is isolated from the lateral sample boundary by 
an annular vortex-free superflow and confined within a 
radius Ry = \J kN/2ttQ,. These central sections of the 
sheets amount to a total length of nRl/b. At the sample 
boundary the 2M ends of the M vortex sheets cross the 
vortex-free region of width d = R — Ry R. We have 

L = ttRI /h + 2Md K 7ri?V6 + 2d [M ~ A/,nax) ■ (3) 

In Fig. 1^, loop (a) is traced by a single sheet (Af = 1): 
(i) In the lower branch of the loop which starts from 
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FIG. 4. Dynamic responses of two different vortex 
sheet configurations to a step change in rotation velocity 
AO, = -0.15 rad/s (at t = 0), with Q = 1.25 rad/s ini- 
tially. The satellite-peak heights (solid lines) normalized to 
the initial values are recorded as functions of time t. The 
thin broken lines show exponential fits to determine the re- 
sponse time T. The thick dashed lines show the results from 
numerical simulations {a/ps = 1.1 • 10~*cm^s~^, qo = 30). 
The initial increase in the uppermost curve is the out-of-phase 
response of the single sheet (at A'^ = const), which becomes 
experimentally visible only at higher drives fl > 1.8rad/s |J13|. 



1.1 rad/s in increasing Q, the counterflow velocity at the 
boundary is at the critical value and new vorticity is cre- 
ated while d is almost constant. We find from Eq. (||) 
that L « TrRl/b cx A^fi"^/^, if we neglect the small 
term 2Md. Since N x fl, the length is i oc fl'^/'^ and 
dL/d^ K, (2/3)L/r2 > 0, which agrees with the posi- 
tive slope above 1.1 rad/s. (ii) During decreasing SI, N 
first remains constant, before the annihilation thresh- 
old is reached. In this branch L oc Sl"^/"^ so that 
dL/dfl « — (1/3)I//J7 < 0. This again agrees with the 
almost twice smaller negative slope when SI first de- 
creases below 1.2 rad/s. The out-of-phase response {L 
increases while the drive S) decreases) is the most dis- 
tinct characteristic of the single-sheet response at con- 
stant number of circulation quanta (Sli < Sl^) fl^ . (iii) 
At the annihilation threshold (branch with S7 < 1.1 rad/s, 
dQ/dt < 0), vorticity is expelled, d w const, and again 
dL/dVl fa (2/3)L/ri > 0. (iv) Finally, the negative slope 
of dL/dfl is resumed in the last leg of the loop where S7 
is increased back towards the critical velocity. 

In contrast multiple sheets. Fig. ^ loop (b), with M sa 
-M^max, have always L oc b'^, as seen from Eq. (||), and 
dL/d^ « (2/3)L/S7 > in agreement with the measured 
narrow hysteresis loop. The in-phase response dominates 
the output from states with multiple sheets. The promi- 
nent difference in the responses of the two sheet configu- 
rations in Fig. ^ has thus a simple explanation in terms 
of one adjustable variable A/. 

Dynamic response - To illustrate the faster dynamics 
of a state with multiple sheets, Fig. ^ provides a com- 
parison to the single sheet. The responses are recorded 
to a step-like change in the rotation drive. The vortic- 



ity relaxes with time constant r = {2aVl)~^ which is an 
order of magnitude smaller for the multiple-sheet state 
than for the single sheet. 

Simulation - These and other measurements can be 
compared to numerical simulation, by extending the 
model in Rcf. jl^ to dynamical processes with an 
anisotropic mobility p^ . In comparing the responses 
of different configurations to slow changes in S7 (Fig. |^), 
it appears that the most important factor is the number 
of sheets, whereas the detailed arrangement of the sheets 
is less influential. Our dynamical simulations confirm 
that in an oscillating drive, a single sheet transformes 
into multiple sheets: When f2 increases rapidly azimuthal 
sections of the sheet become unstable. When the num- 
ber of circulation quanta grows, an azimuthal outer sheet 
becomes overfilled as the quanta have no time to move 
inside the sheet. This section then becomes corrugated 
and a bulge is formed, to redistribute the vorticity more 
evenly. During the subsequent decrease of S7 the vorticity 
expands and the bulge may touch the sample boundary. 
This produces two new connections and improves the ra- 
dial alignment of the sheets. Ultimately, the faster mo- 
tion of vorticity in such sheets suppresses the formation 
of new bulges when M has grown to (0.6 0.8)Mt„ax- 
In Fig. ^ the step response of such a dynamically created 
sheet configuration is compared to that of the single sheet 
from which it was created. In both cases the measured 
curves are in reasonable agreement with the calculations. 

Conclusion - Superfluid "^He-A provides the first exam- 
ple of a quantum system where the temporal properties 
of the external field determine in what form of quantized 
vorticity the response occurs. As a function of frequency, 
a topological transition takes place from linear to pla- 
nar order-parameter configurations. Among all possible 
states of vorticity in •^He-A, a state with radially ori- 
ented multiple vortex sheets is established that provides 
the fastest dynamic response at high frequencies. Many 
of the features, which make this phenomenon possible are 
also present in other systems with a multi-component or- 
der parameter, and thus this could be a generic property 
of quantum systems with intrinsic anisotropy. 
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